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1.        Introduction  a 

The  objectives  of  this   study  are  to  develop  a  numerical 
method  for   estimating  tidal  constituents  in  the  open  ocean  using 
fundamental  theoretical  concepts g    to  adapt  the  method  to   com- 
puter use,    and  to  apply  it  to   the  analysis  of   the  tidal  regime 
in  the  Gulf  of  Mexico « 

The  ability   to   predict   tidal  behavior  in  the  open  sea  is 
useful  in  oceanographic  research  and  coastal  engineering.      It 
permits  removal  of  tidal  effects  from  deep-sea  observations 
and  prediction  of  the  effect  upon  the  tides  of  alterations  of 
the  coast  line0 

Present  tidal  knowledge  is  largely  limited  to  the  results 
of  harmonic  analysis  at  selected  reference  stations  and  empirical 
lag  and  amplitude   corrections  for  closely  adjacent  stations.    This 
technique  is,    of   course,    limited  to   the  coast,    but  the  work  of 
Dietrich  [l]    r  subjectively  fitting  a  two-constituent   (M2  &  Kl) 
model    (see  below)    to  Island  and  coastal  date  between  coastlines 
extends   the  results  into   the  open  ocean.      The  accuracy  of 
Dietrich's  work  is   considered  good  In  dense  data  regions,    such 
as  the  North  Atlantic g    but  In  sparse  data  areas  no  general 
agreement  has  been  reached  concerning  tidal  regimes.      The  dis- 
tortion of  tides  In  shallow  water   cannot  be  accurately  assessed. 
Thus   the  results  of   the  method  are,    to   some  degree,    in  error. 

A  more  satisfying  approach  is  to  apply  the  fundamental 
hydrodynamic  equations  and   calculated  luni-solar  forcing  function 


to  as  realistic  an  ocean  basin  as  possible. 

Two   methods   of  attack  are  the  integration  of   the  basic 
equations  and  the  approximation  of  the  equations  by  finite  dif- 
ferences,   leading  to  a  solution  by  numerical  methods.      In  either 
case,    a  boundary  value  problem  results  in  which  increased  real- 
ism is  accompanied  by  increased  mathematical  difficulty,    where- 
as   oversimplification     can  render   the  problem  physically  meaning- 
less.     In  this  paper  discussion  will  be  confined  to  numerical 
approaches.      Hansen  £2)   used  finite  differences  in  his   treatment 
of  the  M2   tidal  constituent  in  the  North  Sea  with   excellent 
results,    both  in  tidal  heights  and  currents,    as   shown  in  Figure 
1.      Friction  was  included  and  the    beta-plane  approximation  was 
used.      Rossiter's  approach  [3]     involves   the  use  of   spherical 
coordinates  and    frictionless     flow  in  a   meridionally  bounded 
ocean  of   constant  depth.       Coriolis     effects  necessitate     spec- 
ial care  near  the  poles  and  the  parameters  were   modified  to 
insure  convergence  in  these  areas.      In  another  application  lat- 
itudinal boundaries  were  imposed  to  approximate  the  Atlantic 
Ocean.      Relaxation  techniques   were   employed  successfully  by 
both  Hansen  and  Rossiter.      Their  results  and  Rossiter's  recom- 
mendation that  the  method  be  applied  to  irregular   shaped  basins 
were  important   motivations  for  the  present  study. 
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2.        Tides  in  the  Gulf  of  Mexico 

The  tides  in  the  Gulf  of  Mexico  are  modest  in  range,    but 
diverse  in  character  „      Along  the  Gulf   Coast  the  tides   change 
abruptly  and  assume  many  forms.      Figure  2  illustrates  the 
distribution  of  these  regimes «      Along  the  entire  Atlantic   coast 
of  the  United  States  tides  vary  widely  in  range,   but  there  is 
no  variation  in  semi-diurnal  character . 

In  the  light  of   this  constancy  the  abrupt  changes  in  the 
Gulf  are  very  dramatic .      One  example  of   considerable  variation 
of  the  tidal  regime  is  in  the  northeast  portion  of  the  Gulf. 
At  Cedar   Keys,    a   mixed,    predominantly  semi-diurnal  regime 
exists.      The  tidal  curves  vary  markedly  in  both  directions  along 
the  coast.      At   Pensacola,    200   miles  west  of  Cedar   Keys,    the 
diurnal  influences  dominate.      Similar  changes  are  encountered 
in  other  areas  along  the  Gulf  Coast.      A  comprehensive  descrip- 
tion of  the  tides  along  the  Gulf   Coast  and  the  criteria  for 
regime  classification  have  been  given  by  Marmer  (jjj    .      Many 
explanations  have  been  advanced  for  these  shifts  in  the  basic 
form  of  the  tides .      Harris  postulated  in  1897   that  the  semi- 
diurnal tides  were  the  extension  of  the  Atlantic  regime,    while 
the  diurnal  tides  resulted  from  the  natural  period  of   the  Gulf 
basin  o 

Grace  [5]    compiled  the  various  theories  and  applied  an 
inventive  analytical  technique  in  1932.      He  divided  the  Gulf  into 
five  sections  in  each  of  which  the  depth  was  treated  as  a 
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constant.      The    linearized,    inviscid  hydrodynamic   equations   of 
motion  were  solved  within  each  section.      Continuity  was   satis- 
fied at  the  interior  boundaries  between  sections  and  the  bound- 
ary conditions  of  no  normal  flux  were  met  at  the  coasts.    The 
studies  resulted  in  a  unique  set  of   corange  and  co tidal  lines. 
Based  on  this   work  Grace  put  forth  a  theory  of   tidal  dynamics 
In  the  Gulfs    the  semidiurnal  tidal  wave  enters   the  Straits   of 
Florida,    moves  around  the  Gulf   coast  in  a   contra  solem  amphi- 
dromy  and  exits   through  the  Yucatan  Strait  six  hours  later  J 
the  diurnal  tidal  wave  enters  the  Straits  of   Florida  and  prop- 
agates in  such  a  way  as  to  arrive  at  most  points  along  the 
coast  nearly  simultaneously. 

Several  topographical  features  of  the  Gulf   may  well  have 
discernable  effects  on  the  tides.      The  continental  slope  off 
the  west  coast  of  Florida  is  very  steep  and  could  act  as  a 
reflective  surface  to  an  incoming  tidal  wave.      The  intrusion 
of  the  Desoto   Canyon  toward  the  northern  coast  and  the 
Mexican  Basin  are  both  relative  deeps   with  unknown   effects  on 
the  natural  period  of  the  Gulf.      The  Florida   Keys  and  the 
broad  shelf  off  the  west  coast  of  Florida  could  mask  influences 
of  tidal  waves   entering  the  Straits  of  Florida. 
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3.         The  Theoretical   Development 

Introduction 
To   describe  the   tides   utilizing   the  basic   hydrodynamic 
equations   one   must  resort   to   describing   the   separate  tidal 
constituents.      This   limits    each   solution   to  a   particular  known 
period  and  provides   for  a   complete   solution  through  a   summation 
of   constituents.       Except   in  specialized   cases    this   summation 
limits    the  problem    to   linearized   equations   of    motion.    The   work 
described  below   is   directed  towards  a   solution  for   the  principal 
lunar   semidiurnal   constituent,    M2.      Coastal  observations  indicate 
that  knowledge  of   this    constituent   should  permit   explanation  of 
these  observations   and  description  of   regimes   in   the   open  sea. 
Harmonic  analyses   presented  by  Marmer  indicate  that  M2   and 
Kl,    the  lunisolar   diurnal   constituent,    are  the  dominant    constit- 
uents.     The  variability   of    the  amplitude  of   Kl   is   from    .59 
feet   to    .09   feet,    ignoring   the  results   of   short-term   observa- 
tions.     The  amplitude   of   M2   varies   from    zero    to    1.07   feet  by 
the  same  criterion.      Thus    the   spacial  distribution  of   the  M2 
constituent   should  be  indicative  of   the   distribution  of   the   forms 

of   the  tides. 

The   Equation   of   Motion 

From   the   equations   of    motion  and   continuity  we   may  derive 

a   second-order    elliptic   partial   differential   equation  relating   tidal 

height,    depth,    and    their   derivatives.    The  analysis   is   performed 

on  a  beta-plane.    Advective  and  frictional   terms   in  the   equation 

of   motion  are   neglected.    The  hydrostatic  approximation  is   made. 
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The  beta-plane  assumption  restricts   the   size  of   the 
basin  in  which   the  results   will  apply,    but  is   considered  accep- 
table in  the  Gulf. 

Friction  is  present  and  does  influence  tidal  distributions, 

especially  in  shallow   water.    Frictional   effects  are,    however, 
believed  to  be  highly  non-linear,    thus  further   complicating  the 
equations.      (The  work  of   Dronkers£6]  in  linearizing  this   term 
for   shallow   water  is  an  example  of   these  complications.)    Fur- 
ther justification  for  neglecting  friction  is   that   only  the  per- 
imeter is  in  shallow   water  while  most  of   the  Interior  is  in 
relatively  deep  water. 

With  the   stated  assumptions   the  dynamic   equations  are: 

s^-tir^         (l) 

where, 

u  Is  the  velocity  component  In  the  x  direction   (positive 

East) , 

v  Is  the  velocity  component  In  the  y  direction   (positive 

North), 

F  is  the  Coriolls  parameter, 

g  is  the  acceleration  of  gravity, 

a     is  the  tidal  height  above  undisturbed  sea  level   (which 
is   mean  low  water)    and  is  positive  upwards,    and 
f    is  the  forcing  function   (see  below). 

The  hydrostatic  approximation  Implies   that  u  and  v  are 
independent  of  depth. 


llf 


The   equation   of    continuity  is; 

!iwnjLM+if,o  (2) 

where  h  is   the  basin  depth. 

Since   the  frequency   of   the  desired  harmonic   constituent 
is   known,    sinusoidal  forms   for   the  dependent   variables  are 
assumed: 

r-    8e[f(*tA  ei(rt] 

■where  the  notation,    Re  j     ^     indicates   the  real  part  of   the 
quantity  in  brackets, 

U       is   the  angular   frequency   of   the   constituent  under    exam- 
ination,   and 

H   is   the  amplitude   of   the  forcing  function. 
Solving  the  equations   of   motion  for   u  and  v  and  substitut- 
ing the  results   into   the  continuity   equation  yields   one  partial 
differential   equation  for  the   tidal  height.      Writing    ]°  -   f,  *  x  J«- 
( where,    for  instance,//    is   the   component  in  phase  with   the 
forcing  function)    and   separating  the   differential   equation  into 
real  and  imaginary  parts   gives 

-  ^z  /  r  u)  +  ib  $JJ^  -  £  aji  4£  +  ifc  %&'") 

.  £.  ah.  a  A  +  o^l4-  ^  -  o 
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where  ^y      ■= 


=  11  +a* 
a**    s^p  • 


Rossiter  used   spherical   coordinates  in   treating  basins 
covering  large  areas   of  the  earth.      Thus  his   equations  differ 
in  detail  at   this   point  although  the  form  is  very  similar. 

Hansen's   work  in  the  North  Sea  included  friction.      The 
added  terms  in  the  equations  of   motion  had  the  form, 

f%--  R  u       and     FY  =   R  V  (5) 

where  R  is  a  suitably   chosen  proportionality  factor.      The 
excellence  of  his  results  is  a  strong  argument  for  this   some- 
what ad  hoc  hypothesis.      In  the  situation  involving  considerably 
greater  bottom  relief   this  form  is  inadequate  due  to   the  var- 
iability of   the  friction  effect  with  depth.      In  deep  water 
friction  will  not  influence  tidal  distributions  as   much  as  bottom 
topography  which  enters   the  equations   through  continuity   con- 
siderations . 

Boundary  Conditions 
The  boundary  conditions  determine  the  velocities   normal 
to   the  boundaries.      At   solid  boundaries   these  velocities  must 
be  zero.      At  the  passages  into   the  Gulf  values   must  be  assign- 
ed to   the  associated  tidal  currents.      Following  the  treatment  of 
tidal  heights,    we  write 

U   =    UKI     +  I  UK,2 

(o) 

V  -   VKI     +  1  VK2-. 
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Substituting  these  relations   into   the  solutions   for   u  and  v  fi*om 
the   equations   of   motion  and   separating  into   real  and  imaginary 
parts  gives 


^)    [' 


VK2*  life^S^  ^   J 


(7) 


Observations   of   tidal   currents   to   determine   their   con- 
stituents are  not  frequently   made.      For  accuracy  continuous 
observations  are  required   over  long  periods   of   time.    Tidal 
current  measurements   made  by  Harris  in  1887  and  supported 
by  Grace's  work  were  used  in  the    Yucatan  and  Florida   Straits 
in  this  study o      The  accuracy  of  these  calculations  is  question- 
able but  the  order   of  magnitude  is   correct.      By  allowing  flow 
through  these  channels   the   effects   of  a   co-oscillating  tide  are 
included  in  the   solution . 

The  Forcing   Function 

The  forcing  function  is  basically  the  gravitational   effect 
of  the  ideal  body  associated  with  the  constituent  to  be  examin- 
ed o      It  is   obtained  from    equilibrium   tidal   theory    (for   example 
see  ProudmanQ7]  )    in  the  form, 

f°£    M^Lpcos*   peas*   cfcs>szf{( 
H  -   <>>  (  d  )     Me      l* 
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for  the  lunar  semidiurnal  constituents,    and 

H^CL(a)     Ml  2.4  ->       (9) 

for  the  lunar  diurnal  constituents. 
Here, 

a  is   the  radius   of   the    Earth 

d  is   the  distance  from   the  Earth  to   the    Moon 
Mm  is   the  mass   of   the  Moon 
Me  is   the  mass   of   the  Earth 

y     is   the  latitude  of   the  point  under   consideration 
Q       is   the  declination  of   the  ideal  body 
T      is   the  longitude  relative  to   the  reference  longitude. 
For   computations   the   central   meridian  of   the  Gulf,    90°W, 
will  be  used  as  a  reference. 

Similar   expressions  apply  for  solar  diurnal  and    semidiurnal 
constituents  and  also  for  those  with  periods   of  up  to   two   weeks 
In  this   work  these  constituents  have  relatively   small  amplitudes 
and  so  are  not   treated. 

Over  the  area  of  the  Gulf   the  lunar   semidiurnal  forcing 
function  varies  by  about   eighteen  percent. 
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4 „         Numerical   Techniques 

Finite   Difference  Approximations 
To   put   equations    (4)    into  a  form    more  suitable  for   num- 
erical  treatment   the  derivatives  are  approximated  by   central 
differences, 

/^z\  -  San  -"£*->       (d*\  -  ?v-n :  z*-i 

1  51  Jo~      26  ,     igTy  /o"  2€  (10) 


and 


fv^)o  *  I*  ^2xfi+^-i^^t^r.»-4lcj   (11 


) 

whereZ  is  a  dummy  variable  and  the   subscripts   refer   to  loca- 
tion in  the  grid  arrangement  as   shown  below9. 

r 

w-i  • — ±- 


*-£  — «■ 

7T^ »*+l 


Central  differences   utilize  values    sampled  at  an  interval   of   two 
grid  distances o      This  has  a   smoothing   effect  on  the  data. 
Because  the  bathymetry  is  normally  given  to  the  nearest  fathom 
smoothing  is  beneficial  in  decreasing   errors  in  depth  gradient 
due  to   sounding  inaccuracies ,      The  tidal  heights  are  much  smaller 
than  the  depths.      Because  of   this?    solutions  are  very   sensitive  to 
sounding   errors.      For   this   reason  central   differences  are  used 
in  all  finite  difference  approximations. 

The  use  of   central   differences  at  a  point   on  the  bound- 
ary introduces  values   for   tidal  heights   outside  the  boundary  into 
the   equations  o 

For  an  interior   point,    equations    (4)    with  finite  differences 

become 
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and  <12'  I 

(it)      (jlho  +  nx+i-W*-.]  +  (J?)^  (W-  k*+<  +^x-'N)  1 

+  (Jt)^,   C4lno+hHv,-h4.,>)    ^(Jt)y.i  (4^0-hvn+Kx-^  I 

^(^eC-^M*  (*T-n)wE£(K+>-^3  ! 

Consideration  was  given  to  variable  grid  distancing^  to  permit 
detailed  analysis  in  areas  of  large  depth  gradients.      However, 
the  stability  of  the  computation  about  a  point  with  different 
distances  to   each  of   the  adjacent  points  is  questionable,    and 
the  finite  difference  form   of   the  equations  is  quite  complex. 
Appendix  B    contains   the  finite  difference  equations  for  a 
variable  grid  distance.      The  density  of  bathymetric  data  in 
the  Gulf  is  such  that  a  grid  finer  than  that  used  below  would 
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require   considerable  interpolation  and   thus   no   improvement  in 
the  solution  would  appear. 

Rossiter   used  a  finite  difference   equation  involving  der- 
ivatives at   the  boundary  and   three  interior   points   to   produce 
the  tidal  heights   outside  the  boundary.      Although   it  is  valid, 
Rossiter's   procedure  has   the  disadvantage  of   introducing  para- 
meters associated  "with  non-adjacent  points  into  the  equations. 
This   complicates   the  system   of    equations   in  a  relaxation  tech- 
nique by  increasing   the  number   of   points  involved  in  each   equa- 
tion by   two.      It  also   disrupts   the  iterative  form   designed  for 
computer  use,    which   will  be  discussed  later. 

At  a  boundary  point,    the  basic   equations    (12),    derived 
from   continuity  considerations   and  the  applicable  boundary   con- 
ditions are  both  in  effect.      These  boundary   conditions    manifest 
themselves   in  the  forms   of   prescribed  velocities   normal   to   the 
boundaries. 

The  velocity   equations    (7)in  finite   difference  form  are 

(13) 
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Thus,    for  a  given  normal    velocity,  the  tidal  height  for   the 
point  outside  the  boundary   may  be  solved  for  in  terms  of   three 
interior   tidal  heights.      Substituting  these  solutions  into   the 
equation  for  an  interior  point  now  incorporates   the  dynamic 
equations,    the  continuity  equation  and  the  boundary  condition 
into  a  set  of   two   equations  at  the  point,   as  well  as   eliminat- 
ing values  outside  the  boundary  from  the  calculations. 

The  substitutions  and  resultant  equations  will  be  unique 
for  each  boundary   or   corner.      As  an  example,    to   solve  for  a 
point  on  the  south  boundary  with  a  known  velocity,    equations 
(13)    are  used  in  the  form, 

fjr-/vjY„  -  rr-«)Y+-|(/a)x.,^(jiv.,-v^^) 

(14) 

and  substituted  into   the  Interior  equations    (12)    to   eliminate 

values  Involving  the   exterior  point  Y-l. 
This  results  in 

+  riT-w)>cMC4Uo-(fh^-h)c-l)fl-FV)J  +  4K0M  <«> 

and 
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To  accommodate  the  boundary   conditions   involved  at  a   corner 
values   of  u  and  v   must   enter   the  basic   equations    ( 12 )  .    This 
is  accomplished  by   simultaneous   solution  of   equations    (13)    to 
eliminate  terms   involving  the  two   heights  at  the  same  point 
from   the  combined   equations.      The  results  are 

(16) 


Combinations   of   these   equations   will   meet  any  corner  boundary 
condition.      For   example,    substituting  to   eliminate  values  at  the 
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point  X+l   in  the   south  boundary    equations   will  give   equations 
for   the   southeast   corner   in  the  form, 

--  o 


(17) 


and 


The  forms   of   the   equations   for    each  boundary  and   corner 
point  are  presented  in  Appendix  A. 
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Methods    of    Solution 

When  equations    (12)    and  their   equivalents   incorporating 
boundary   conditions    (see  Appendix  A)    are  applied  at   each  point 
of  a  grid   consisting  of   n  points   describing  the  body  of   water 
studied,    a   system   of  2n   equations  in  2n  unknown     tidal  heights 
results.      Both   methods   of   solution  undertaken  here  are  pred- 
icated upon  the  presence   of  a   dominant  diagonal  in   the  system's 
coefficient   matrix.      That   is,    in  the  application  of   the  equations 
to    each  point  there  is   one   coefficient   much  larger   than  all 
others  in  the  equation  and  in  each  successive   equation  the 
largest   coefficient  multiplies  a  different  unknown  value  of  the 
tidal  height.      Both  methods  are  designed  to  obtain  a   solution 
by  altering  the  unknown  tidal  height  associated  with   the  largest 
coefficient. 

Gauss   -    Seidel  iteration   ( see  for    example,    Hildebrand  LAJ   ) 
is  the  method  used  in  the  computer  program.      Briefly,    this 
method   consists   of   obtaining  new     values  for   the  variable  with 
the  dominant   coefficient,    say  X^    ,   In  an  equation  of   the  form, 

A,  X,   +Aatf*+   •  •  •+  AvmXm  "-  0        («) 

by  the  formula 

This  satisfies  the  equation  with  approximate  values  of  X.-,...., 
Xn  by  varying  X^ only.  This  procedure  is  carried  out  for  each 
equation  in  the  system   following  a  first  guess   of   the  value  of 
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< 

each  variable.      Repeated  use  of   equation   (18)    usually  brings 
about   "convergence",    in  the  sense  that   each  of   the  equations 
is   satisfied  by   the  same   set  of   unknowns   to   within  prescribed 
tolerances.      Gauss  -    Seidel  iteration  is  not  highly   sensitive 
to  the  first  guess  values. 

Relaxation  is  a  similar  process  but  permits   selective 
variation  of  the  unknowns   to  bring  about  convergence.    Using 
the  same   example   equation  the  first  guess   will  give 

A.Xi  +  Ai*i+  ••-+   AnXn    -  K 
where  R   is   the   "residual",    which  should  vanish  to  achieve  con- 
vergence of   the  system   of    equations.      Changes   will  be  made 
in  the  unknown  associated  with  the  largest   coefficient  following 
the  rule 

X.'-X.-'N-Sr  (19) 

where  TV  is   the   "Relaxation  Coefficient"  and  may  take  on  values 
in  the  range 

0  <  >  <   z  .  (20) 

When  »^  =1   the  residual  will  be  reduced  to   zero  and  the 
process  is  very  similar  to  the  iteration  described  above.      When 
~h>|    the  residual  tends   to  alternate  sign  on  successive  iterations 
and  decrease  In  magnitude.      This  process  usually  speeds  up 
convergence  and  is   called  overrelaxation.      When~^\the  residual 
tends   to  be  reduced  in  magnitude  and  keep  the  same  sign  thus 
approaching  zero  from  one  side  rather   than  oscillating  around 
zero.      This  Is   called  underrelaxatlon. 
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5«         Hand   Calculation  of   the  M2   Tidal   Constituent  in  the  Gulf 

To  gain  an  understanding   of   the  behavior   of  a   tidal   solu- 
tion in  the  Gulf  and   to   supplement   the  computer   work  a  hand 
solution  was  undertaken  for  a   relatively   coarse  grid.      The  grid 
arrangement   consists   of    16   points   with  a  grid  distance   of   180 
nautical   miles 0      Figure  3    shows   the  arrangement.      This  grid 
represents  the  Gulf  fairly  well,    but  has  the  disadvantage  that, 
in  areas   of  large  depth  gradient,    the  terms  in  the   equations 
may   not   occur   with   the  relative   magnitudes   desired  for  a  rel- 
axation solution. 

Applying   equations    (12)    and  appropriate  boundary   conditions 
to   each   of   the  16   grid  points   produces  a   system   of  32    equations 
for  32  unknown  tidal  heights.      Of  these,    eight  do  not  contain 
the  dominant  coefficients  desired  for  relaxation.      Six  of   these 
eight  are  for  interior  points    with   eight  terms  in   each   equation. 
This  means  that   successful  reduction  of  the  residuals   of   these 
equations,    especially  in  the  later   stages   of  the  solution,    may 
depend  upon  systematic  adjustment  of   the  number   of  tidal  heights 
to  avoid  deleterious   effects   on  the  residuals  at  adjacent  points. 

In  every  iteration  values  of  the  tidal  heights  used  were 
the  results  of  the  last   calculations  made  at  the  point.      Some 
systems    (e0g0,    the  method  of  Richardson  as  described  by 
Hildebrand    [i  bid.]    )    are  cyclic  in  nature  using  only  calculations 
from  previous  iterations.      The  approach  in  this   work  is   not 
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cyclic,    but  is  based  upon  maximum  reduction  of   the  largest 
residuals   consistent  with  the  long-range  effect  upon  adjacent 
points. 

Initial  passes  were  made  based  on  a  first  guess  field  de- 
rived from   Grace's  results  and  observational  data  reported  by 
Marnier.      Over-relaxation  techniques  were  utilized  with  "N  tak- 
ing values  from   1,5   to   1.2.      The  chief   drawback  of  this   tech- 
nique is   that  it  can  produce  oscillations  in  the  values   of   the 
tidal  heights  and  thus  fail  to  converge.      At  the  first  sign  of 
an   undamped  oscillation,  P-  was  lowered  to   1.      On  subsequent 
iterations  the  relaxation  factor  was  varied  to   meet  the  needs 
of   each   equation.      To  determine  corrections  late  in  the  cal- 
culations  systems   of   simultaneous   equations  were   solved  for 
the  appropriate  values  of  1^  .      At  this   stage  the  equations 
were  very  sensitive,    and  large  corrections  frequently  had  del- 
eterious  effects  on  adjacent  points.      To   most   efficiently  deal 
with  these  situations   tabulations  of  present  residuals  and  the 
effects  of  proposed  changes  in  tidal  heights  upon  them  were 
necessary.      These  were   most  helpful  In  determining  the  Inter- 
action of  changes  In  the  tidal  heights.    After  twenty  iterations 
component  tidal  heights   (^Tandj^)   at  all  Interior  and  solid  bound- 
ary points  were  defined  to  an  accuracy  of  plus  or   minus   three 
centimeters.      Those  at   the  Yucatan  Passage  and   Straits  of 
Florida  were  extremely  sensitive  to  values  of  Input  current, 
fluctuating  approximately  ten  centimeters  with  changes   of  only 
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one  centimeter   per   second  in   the  associated  velocity  component. 
The   effects   of   these  points   on  other   equations  is   not  large. 
Since   the  current  velocities  utilized  are  not   considered   suffic- 
iently accurate  to  warrant  it,    these  points  were  not  relaxed 
to   the  degree  applied   to   the  other  points  in  the  grido 

The  results   of   these  calculations  are  contained  in  Appendix 
D,      In  general9    this   solution  was  very  similar  to   that  accom- 
plished using  the   computer   on  the  same   16  point  grid.      The 
maximum   difference  in   magnitude  of   tidal  height   was  9gj6   cms. 
The  average  difference  was  3°23    cms.      This   computer   solution 
is   discussed  in   more  detail  below. 
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6.         Solution  by   Digital   Computer 

Adaptation   of    the  finite  difference  formulae  ('12)    to   the 
computer  is  a   difficult   task.      For    efficient  iteration  a  proced- 
ure is  necessary  for   expressing  all   equations  in  a  common  for- 
mat.     This   procedure  permits   the  use  of  indexing  operations 
in  the  identification  and  location  of   equations.     This   utilization 
of  indexing  is  accomplished  by   the  inclusion   of   identification  of 
adjacent  points  and   their  relative  positions  in  the  input  data. 
Through   this,    and  a   coded   method  for   points   outside  the  bound- 
ary the  program   identifies  interior  and  specific  boundary  loca- 
tions and  simultaneously  applies   the  proper   equations   to   them. 
Fortran  60   program  language  is  used.      Based  on  the  needs  of 
this  particular  application,    convenient  computer   capacity,    and 
the  concept   of   versatility,    the   maximum  number   of   points   in 
the  grid  is   1200.     These  points  are  numbered  serially,    and 
points  outside  the  boundary  are  identified  by  numbers  greater 
than   1200.      Since   the  requirement  for   ficticious   tidal  heights 
outside  the  boundary  has  been  eliminated   (see  above),    the  only 
data  required  outside  the  boundary  are  tidal  current  velocities 
and  basin  depths. 

The  presence  of  an  identification  number  greater   than 
1200   then  specifies  a  particular  boundary  and  two  such  numbers 
specify  a  particular  corner.      The  current  velocities  included  as 
data  for  points  outside  the  boundaries  replace  tidal  heights  for 
those  points  In  the   equations  at  the  boundary  point.      This 
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requires  alteration  of   the  coefficients  in  the   equations   to  fit 
the  location  of   the  point  at   which   the  equations  are  being  ap- 
plied.     This  is  accomplished  from   the  general  forms, 


and 


4-  (Xa  (£-H)  r+i  -  b4  (  Jl )  y-i 

+  cj/  ^ 


(21) 


(fi)o   =  [cuf  /?)*+>  *  b,  (r,-N)^ 

Thus   the  coefficients   depend  upon  the  location  of   the  point. 
Appendix  A  presents   these  coefficients  in  detail  and  the   com- 
puter program   is  given  in  Appendix  C. 

The  system  of   equations  derived  from   the  grid  utilized 
in  the  hand  calculations  was   then  solved  using  the  computer. 
Actual  depths  and  currents  were  employed.      This   system  con- 
verged in  37  iterations. 

In  a  study  of   this  nature  results  are  best  displayed  by 
the  analysis  of  the  solution  to  obtain  cotidal  and  corange  lines. 
Gotidal  lines  are  lines  along  which  the  high  tide  occurs   simul- 
taneously.     Corange  lines   are  lines  along  which   the  tidal  ampli- 
tudes are   equal  .As  set  forth  above  these  results  apply   to   the 
M2   constituent,    and  the  phases  are  relative  to   the  passage  of 
the   moon  through   the  90th  west   meridian.      The  lack  of   detail 


31 


in  this  solution  causes   the  corange  and  cotidal  analysis   to  be 
somewhat  subjective.      The  analysis   of  the  16  point  solution  is 
illustrated  in  Figure  k   which  is  based  upon   coastal   data  as   well 
as  the  computer   solution.      It  is   evident  that  most  of   the 
distortion  of  the  tides  occurs  in  the  coastal  regions.      This 
does  not   mean  that  the  solution  is  In  error.      The  lack  of 
detail  does  not  permit  identification  of  a   small-scale  amphidromy 
such  as  that  proposed  by  Grace  or   seen  In  the  solution  discussed 
below.      The  path  of  the  progressive  tidal  wave  Indicated  by  the 
cotidal  lines  is  logical  and  amplitude  distributions  agree  with 
observed  coastal  data. 

To  obtain   more  detail  In  the  solution  a  finer  grid  was 
desired.      Using  a  grid  distance  of   thirty  nautical  miles  a  grid 
was   then  constructed  In  an  Irregular   shape  closely  approximating 
the  Gulf.      Real  depths  were  used,    taken  from  U.    S.    Coast 
and  Geodetic   Survey  Chart   1007.      The  grid  consisted  of   504 
points.      Current  data   mentioned  above  were  used  in  the  Straits 
of  Florida  and  Yucatan  Strait.      The  resulting   system   of   equa- 
tions failed  to  converge  after  2500   iterations.      In  the  initial 
iterations   the  tidal  heights  increased  slowly  In  magnitude. 
After  25  iterations   the   solutions  were  unrealistic  and  they 
Increased  more  rapidly  In  magnitude.      Computational  Instability 
seemed  to  be  present  and  further   efforts  were  directed  at  a 
more  tractable  situation. 

Using  the  same  grid  a  constant  depth  of   1000   fathoms 
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was   introduced.      The   system   of    equations   converged  in  378 
iterations   with  a  relative   convergence   criterion  of  0.1.      That 
is,    each  tidal  height   within  the  grid   must   change  by  less   than 
10$   of  that   of   the  previous   iteration  before   the  iterative 
procedure  is   stopped.      The  maximum   change  of   tidal  height 
over   the  last   three  iterations   was    .12   centimeters.      These 
results  are   tabulated  in  Appendix   D.      Figure   5   shows   the 
analysis   of   these  results   in   terms   of   corange  and   cotidal  lines. 
The  path  of   the  M2   constituent  tidal  wave  crest   through  the 
Gulf  is   described  by  the   cotidal  lines  labeled  in  hours  with 
zero  being  the  time  of   the  passage  of   the   moon  through   the 
90th  west   meridian. 

The  detail   obtained  in  a   504  point  grid  considerably  lessens 
subjectivity  in  these  cotidal  and  corange  analyses.    Significant 
aspects  of  these  analyses  ares    the  well-developed  amphidromy 
in  the  Northeastern   corner   of   the   Gulf,    the  presence  of   the 
greater  amplitudes  around   the  coast,    the  presence  of  a  rather 
flat  field  of  low   amplitudes   in  the  vicinity   of   the  amphidromic 
point,    and  the  relatively  high  amplitudes  in  areas   of  known   semi, 
diurnal  dominance.      The  path   of   the  tidal  wave  crest   through 
the  Gulf  is   very  well   defined. 

Both  the   16-  and   502+-  point  grid   computer   solutions 
produce   corange  and   cotidal  lines   consistent  with   coastal   data. 
Because   one  is  based  upon   constant   depth  and  the   other  upon 
actual   depths    the  analyses   differ   considerably  in   some  areas ; 
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however,    the  regime  shifts   mentioned  in  Section  3    (say, 
semidiurnal  to  diurnal)    are  clearly  indicated  in  both  Figures 
4  and  5« 
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7.        Conclusions  and  Acknowledgements 

The  three  solutions   obtained  Indicate  that  the       numerical 
methods  previously  discussed     are  effective  in  predicting  tides 
at  sea.      Based  on  the  relatively  quick  convergence  with  actual 
depths  of  the  16  point   system  which  has  a  grid  distance  of 
180   nautical  miles,    the  convergence  with  constant  depth  of 
the  504  point  system   which  has  a  grid  distance  of  30   miles, 
and  the  failure  to  converge  with  actual  depths  of  the  504 
point  system,   it  is  reasonable  to  assume  that  there  exists  an 
intermediate  grid  distance  such  that  convergence  may  be  ob- 
tained with  maximum  detail  in  the  solution.      Another  approach 
is  to  obtain  the  maximum  detail  within  the  field  by  mathemati- 
cally smoothing  the  actual  depth  data  to   such  a  degree  that  the 
504  point  system  -would  converge.      The  danger  here  is   that 
excessive  smoothing  could  distort  the  data  thus  giving  less 
realism  in  the  solution. 

The  cotidal  and  corange  analyses  of  the  constant  and 
variable— depth  models  Indicate  that  variations  in  depth  are  not 
solely  responsible  for  the  tidal  regime  shifts  briefly  described 
in  Section  1,  (see  Figure  4).  Additional  studies  may  deliniate 
the  relative  Importances  of  spaclal  forcing  function  variation, 
input  current,  basin  shape,  and  basin  depths  upon  tidal  ampli- 
tude distributions  and  perhaps  uncover  other  important  variables. 

By  removing  certain  topographical  features   of   the  Gulf  floor 
from   the  data  It  would  be  possible  to  assess  their   effect  upon 
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tidal  distributions.      Similarly  variations  in  input   currents  and 
the  forcing  function  field  could  be  examined  for   their   effects. 
By  adjusting  input  currents   only,    under  a  condition  of   maximum 
realism  and  detail  in  all  other  data,    to   obtain  the  constituent 
distribution  most  compatible  with  coastal  observations   one  would 
be  able  to  evaluate  the  accuracy  of  present  constituent  tidal 
current  data. 

The  only  geographical  limitation  in  this  technique  is   the 
size  of  the  basin  due  to  the    beta-plane  assumption.      Thus  all 
other  seas  comparable  in  size  to  the  Gulf  of  Mexico  are  well 
suited  for  an  investigation  of   this  sort. 

This  study  and  the  writer  have  both  benefitted  greatly 
from  the  Interest,    guidance  and  considerable  abilities  of   mem- 
bers of   the  faculty  and  staff  of  the  United  States  Naval 
Postgraduate  School.      Principal  among  these  are:   Associate 
Professor  T.    Green,    III  for  his   encouragement,    help,    con- 
siderable effort  and  constructive  guidance}    Professor  G.    «J. 
Haltiner  for  his  counsel  and  personal  example  of   scientific 
excellence;   Associate  Professor  J,   B.   WIckham  for  his  Incisive 
and  constructive  criticism}    R,    D.    Brunell     for  his  personal 
Interest,    efforts  in  the  programming  and  technical  assistance; 
also   Dr.    G.    B.   Austin  and   Dr.    G.    Dowling  of   the  U.    S.    Navy 
Mine  Defense  Laboratory,    Panama  City,    Florida  for  their 
interest  and  assistance  in  selection  of  this  topic;   and  Mrs. 
«Julene  D.   Gainer  for  her  patience,    encouragement,    understand- 
ing and  recently  acquired  proficiency  at  the  desk  calculator. 
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Fig.  i'.  Representation  of  amplitude  and  phase  of  the  M2  constituent  in  the  southern 
North  Sea.  The  length  of  the  arrows  is  proportional  to  the  amplitude,  and  the  direc- 
tion is  determined  by  the  phase.  Full  arrows  after  computation ;  broken  arrows  after 
observations.       From  Jf*nt  «ft    {2)  * 
26- s.  1. 


1 


38 


Figure   2    -   GULF    OF    MEXICO    TIDAL    REGIMES 
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Figure  3.    -   GRID  ARRANGEMENT    FOR    HAND 
CALCULATIONS 

This  consists  of   16  pc  jrid  with   180   nautical 

mil©  grid  distance.      Input  currents  were  u   at     point  Ik 
and  v  at  point  13. 
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Figure   4.    - 


RESULTS    OF    VARIABLE    DEPTH 
COMPUTER    SOLUTION: 


Corange  lines    (dashed)    and    cotidal   lines    (solid). 
Corange  lines   are  labeled   in   centimeters.    Cotidal   lines 
are  labeled  in   degrees   denoting   the   phase   of   the   tide 
relative   to    meridional  passage   of   the   moon   through 
90°W. 
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Figure   5.    -  RESULTS    OF    CONSTANT    DEPTH 
COMPUTER    SOLUTION: 

Corange  lines    (dotted)    labeled  in  centimeters, 
Cotidal  lines    (solid)    labeiec     n   hours   with   zero  being   the 
time  of   the  moon's  passage  through  the  90th  West 
m  eridian  • 
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APPENDIX  A 

FINITE  DIFFERENCE  EQUATIONS  FOR  INTERIOR 
AND  ALE  BOUNDARY  POINTS  ADAPTED  FOR 
COMPUTER  SOLUTION 

Equation    (12)    represents   the  finite  difference  form   for 
an  interior  point.      We  have   shown  that  at   each  boundary  un- 
ique  modifications   of   this    equation  are  necessary  to  represent 
continuity    and  the  boundary   conditions.      To  avoid  the  necess- 
ity of   separately   establishing   each   equation  within  the  program 
and   to   create  a   compact  iteration   scheme  the  following   system 
was   devised. 

Examination  of  all   the   equations  led  to  a   general  form, 

+  dz(fr-«)^     +  b3(7l)Y+i 

-r    C    /   C2. 

where  the   symbol  /     indicates   that  all   terms   on   the  right   side 
are  divided  by  C2 

and  f  ^       \ 
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These  equations  are  set  up  to  carry  out  Gauss-Seidel 
iteration,  but  modification  for  other  solution  methods  could 
be  made  at  this  point. 

Now  by  indexing  within  the  program   the  constants   may 
be  called  to  fit  any  of  the  previously  derived   equations   satis- 
fying both  dynamic   considerations  and  boundary  conditions.    For 
interior  points   the  constants  are: 

For  the  West  Boundary, 

Cr   K4'  V  Ci--i(oh.-  4£lCa^F_^ 

East  Boundary, 

Gt5  Bho  bts    O 

cu  ~  4Wo+(hiA+rh^J(l-%y     b5-  -  "5^" 

4f*  frr»-F^  H        /V  '  t/   ^-4^^^ 
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South  Boundary, 


North  Boundary, 


-4hoF 


at^ino-CU^-Kx-»)0"  /cr  )    b*-      o* 
04"-  8ho 


2e  ftr'-F')  ('4-U»<-h4ti-li^-'') 


8 


Southwest  corner, 


Ai 


ft  h  cfr»-f  »^    y 


a 


bi   -     o         ,     K    , 

bi   -      ocrO-FV^  ^  ' 

b^   *     O 

Cz  »   IG  ho-  *  " 


Northwest  corner  ,v 

Or-  Sho 

C    -    46xH«rU»j 


b. 


-2€fox-?2\ 


V 


(tk*^-^) 


bi  * 

b4  *    o 


Southeast   corner, 


C\  L 


;  s ,  bz  -   o 


Northeast   corner, 


Ci-  ^  Ci-  l^ko* 


*& 
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APPENDIX    B 


DIFFERENCING    EQUATIONS    FOR   VARIABLE 
GRID    DISTANCE 

Given  a   relative   grid  arrangement   of    the   sort, 


AY< 


X-l     - 


^u 


-AX* 


%+\ 


AY. 

it-   is   desirable   to   use  only   contiguous   points   in   the  differencing 
scheme  because  the  use   of   other   points  "would  require  a  field 
of  values  and  points   outside  the  real  grid.      This   in  turn  -would 
negate  the  use  of   central   differences   in   evaluating  the   second 
derivative. 

The  basic   forward  differencing   formulae  were  applied  to 


give 


v2z 


(At-)fe.*+)*-(&)U)(Z*-\ )  -  2     Ifo  <,V ■+  (**S\ 

,  (ay.)  &y+!\  +  (tf+)  (Iy-^)-?_o  Ray,)  +(^-)l 

(AY-)2"  AY+ 


where  2   is  a   dummy  variable.      This   expression  is   quite  involved 
and  the   stability  of   calculation  about   such  a  point  is  highly 
questionable. 
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APPENDIX  D 


THE  NUMERICAL  RESULTS  OF  THE  504  POINT 
CONSTANT  DEPTH  SOLUTION 


POINT 

1. 

2. 

3. 

4- 

5. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 

21. 


MAGNITUDE    (cms)  PHASE    (Degrees) 

1.3887000  -141. 24518 

-138.80641 

-131.90854 

-120.89377 

-105.43202 

-86.843942 

-73.334869 

-I72.507H 

-172.30618 

-171.97713 

-171.56044 

-140.26361 

-137.87719 

-130.61788 

-118.38878 

-99.810982 

-.75.709188 

-53.916750 

-39.475631: 

-23.200350 

-4.1130445 


1.3756521 
1.3217907 
1.2539393 
1.2029531 
1.2059750 
1.2334819 
1.7517582 

1.7194041 

1.7070080 

1.6958950 

1.3027698 

1.2262120 

1.1332772: 

1.0404123 

.98849004 
1.0718649 
1.4355376 
2.4490565 
4.1602358 
7.4439120 
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22. 
23  • 
24. 
25. 
26. 

27. 
28. 
29. 
30. 

31. 
32. 

33. 
34. 
35. 
36. 

37. 
38. 
39. 
40. 
41. 
42. 
43. 
44. 
45. 
46. 


5.4827457 
1.6776942 
1.6361282 

1.5939419 
1.5654488 
I.5564OH 
1.5592573 
1 . 5838993 
1.5963884 
1.5907292 
1.5549091 
1.4619722 
1.2553052 
1.1279148 
1.0039504 
.88153194 
.79735342 
.86428115 
1.2173230 
1.9446390 
2.9908028: 
4.1469681 
3.7697796 
2.9949950 
1.5900268 


-2.1818875 

-172.36824 

-172.15673 

-171.76625: 

-171.39348 

-170.96102 

-170.38945 

-169.61342 

-167.86595 

-165.26372 

-161.59242 

-156.27691 

-148.15227 

-142.15830 

-134.12128 

-131.44064 

-100.55519: 

-71.147404 

-45.294358 

-28.377334 
-15.251722 
-3.9740857 
1.7862483 
9.8396774 
-171.85187 
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47. 
48. 
49. 
50. 
51. 
52. 
53. 
54. 
55. 
56. 
57. 
58. 
59. 
60. 
61. 
62:. 
63. 
64. 
65. 
66. 
67. 
68. 
69. 
70. 
71. 


1.5265542 

1.2+871450: 

1.4817022: 

1.4258745: 

1.2+202776 

1.2+22+2894 

1.2+326960 

1.2+367689 
1.2+264928 
1.39102+75 
1.3153951 
1.1943602+ 
1.0752430 
.94189596 
.79343153 
.64595142 
.59415964 
.80011771 
1.2755715 
1.8856044 
2.4422500 
2.5285018 
2.5386045 
1.5837556 
1.5189584 


-171.55308 
-171.29280 
-171.102+65 
-170.47495 
-169.98678 
-169.33126 
-168.36735 
-166.84673 
-164.70118 
-161.82801 
-158.02929 
-153.24541 
-148.40719 
-142.03557 
-132.02672 
-113.55350 
-79.717909 
-44.336330 
-22.528501 
-9.0438157 
1.3021725 
10.436640 
20.977383 
-171.35588 
-171.12447 
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72. 
73. 
74. 
75. 
76. 

77. 
78. 
79. 
80. 
81. 
82. 
83. 
84. 
85. 
86. 
87. 
88. 
89. 
90. 

91. 
92. 

93. 
94. 
95. 
96. 


1.4459243 
1.3873289 
1.3523482 
1.3203357 
1.3014904 
1.3007594 
1.3081634 

1.3162539 
1.3231527 
1.3207365 
1.3008327 
1.2537422 

1.1791414 
1.0871412 

.9684444 
.82084260 

.63568331 
.42310620 
.32467912 
.55898858 
.93149353 
1.2945355 
1.5414282 

1.8510571 
1 . 4426962 


-170 . 90635 
-170.62673 
-170.37487 
-170.00174 
-169.57599 
-169.12103 
-168.52707 
-167.68322 
-166.52922 
-165.02452 
-163.15918 
-160.91726 

-158.41437 
-155.84092 

-152.64950 
-148.05767 
-139 . 73622 
-119.25154 
-66.502209 
-18.151299 
2.2351259 
15.742440 
29.348534 
43.251027 
-170.27041 
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97. 
98. 

99. 

100. 

101. 

102. 

103. 

104. 

105. 

106. 

107. 

108. 

109. 

110. 

111. 

112. 

113. 

114. 

115. 

116. 

117. 

118. 

119. 

120. 

121. 


1.3801700 
1.3135581 
1.2602636 
1.2306716 
1.2046669 
1.1925103 
1.1990822 
1.2142370 
1.2307824 
1.2495502 
1.2644726 
1.2694225 
1.2565822 
1.2247704 
1.1738689 
1.0947432 
.98741067 
.84121204 
.64327279 
.40957519 
.16316727 
.27612851 
.60679389 
.98720807 
1.4865882 


-170.11656 
-169.87776 
■169.60957 
-169.39582 
-169.06696 
-168.71396 
-168.38716 
-167.98880 
-167.44851 
-166.79710 
-166.04208 
-165.22816 
.I64.40579 
.163.71041 
-163.24897 
-162.92580 
-162.80156 
-162.89192 
.163.50925 
■167.70662 

160.04043 
75.916858 
67.911211 
74.397350 
81. 249956 
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122. 
123. 
124  o 
125. 
126. 
127. 
128. 
129. 
130. 

131. 
132. 

133. 

134. 

135. 

136. 

137. 

138. 

139. 

140. 

141. 

142. 

143. 

144. 

145. 

146. 


1.3192293 
1.2570372 
1.1937553 
1. 1450  593 
1 .  1211817 
1.1020564 
1.0981462 
1.1145  761 
1.1415354 
1.1725985 
1.2102442 
1.2500523 
1.2873929 
1.3153326 
1.3320425 
1.3344956 
1.3135021 
1.2  716266 
1.2020871 
1.0981960 
.9868950 
.88220144 
.86646759 
.98653169 
1.2799964 


-168.98529 

-168.94400 

-168.73472 

-168.49526 

-168.34803 

-168.10229 

-16  7. 88417 

-167.77086 

-167. 67201 

-167 .  54610 

-167. 46968 

-167. 47233 
-167.60798 

-167.92557 

-168. 52434 

-169.47107 

-170.77450 

-172.57444 

-175.05351 

-178.68367 

175.47563 

165.61313 

150.92426 

136.05371 

126.90135 
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147. 
148. 
149. 
150. 

151. 
152. 

153. 

154. 

155. 

156. 

157. 

158. 

159. 

160. 

161. 

162. 

163. 

164. 

165. 

166. 

167. 

168. 

169. 

170. 

171. 


1.7474118 

1.2092843 

1.1461983 

1.0849542 

1.0403817 

1.0224132 

1.0109213 

1.0166413 

1.0451478 

1.0872336 

1.1373998 

1.1995002 

1.2705283 

1.3470828 

1.4228629 

1.49577548 

1.5614789 

i. 6100079 

1.6439044 

1.6564831 

1.6410146 

1.6223832 

1.6032476 

1.6392330 

1.7547253 


121.69064 
-167.64318 
-167. 67550 
-167.48544 
-167. 27501 
-167. 21450 
-167.08640 
-167.05710 
-167.22489 
-167. 49308 
-167.83048 
-168. 32647 
-169.00223 
-169.88820 
-171.00051 
-172.38357 
-174.05059 
-175.98955 
-178.23670 

179.18471 

176.18332 

172.66604 

168.28565 

161.61001 

155.14506 
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172. 
173 . 
174  ■ 
175. 
176, 
177 . 
178o 
179. 
180  o 
181. 

182 . 

183 . 
184. 
185. 
I860 
187. 
188. 
189. 
190  o 
191. 
192 . 
193. 
194  c 

195  o 

i960 


2,0045504 

2 o 3098794 
3.0581112 

1.1086753 
lo  0442983 
o 98475467 
.94421575 
.93236154 
o9290726l 

.94545533 
.98770361 
1.0474125 
1. 1201448 

lo  2111263 

1.3187259 

1.4407064 

1.5713587 

1.7083588 

1. 8462118 

1.9728861 

2o0883663 

2.1825526 

2.2449039 

2.2967387 

2o3377767 


150.74902 
149.75370 
155.39824 
-166.33404 
-166.35663 
-166. 13134 
-165. 92555 
-165. 95960 
' -165.97488 
-166. 17686 
-166.67282 
=»167.33990 
-168.13757 
-169.14457 
"170. 35560 

-171.76804 
-173. 36616 
-175.15459 
-177.11356 
-179.21124 
178,57807 
176.30378 
174.02875 
172,02866 

171.09147 
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197. 
198. 
199. 
200. 
201. 
202. 
203. 
204. 
205. 
206. 
207. 
208. 
209. 
210. 
211. 
212. 
213. 
214. 
215. 
216. 

217. 
218. 
219. 
220. 
221. 


2.4036320 
2.5284900 
2.7636824 
3.1167505 
3.5327726 

1.0129899 
.94788364 
.89049944 
.85424915 
.84869031 
.85389577 
.88147957 
.93842593 
1.0173279 
1.1149818 

1.2381205 
1.3866221 

1.5595657 
1,7520227 
1.9616174 
2.1814858 
2.3965146 
2.6017553 
2.7802197 
2.9135588 


165.81038 

162.15293 

159.89052 

158.85404 

157.76749 

-165.13166 

-164.99387 

-164.63719 

-164.38986 

-164.51702 

-164.69366 

-165.15934 
-166.01340 
-167.08388 

-168.30477 

-169.73308 

-171.33510 

-173.08435 

-174.95054 

-176.92430 

-178.97981 

178.91438 

176.79859 

174.72718 

172.74419 
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223 . 
224. 

&£*■  _^  o 

226o 

228 . 

230  o 

231, 
232 . 

233  o 

234. 

235* 
236. 

237  a 

238  c 

239  c 
240. 

241  „ 

242  o 

243  - 

244  ° 

245  - 

246  • 


3.0161880 
3*0826239 
3 o 1604651 
3 o 27283 19 
3 o 4899034 
3 o84437 13 
4»3693820 

5.232202? 
1. 0182813 
.91611772 
.85362642 
.79976933 
3 76826352 
.76900960 
c 78258238 
o 82125881 
.89298309 
.99155335 
1.1152134 
1.2724424 
1.4649441 
1 o6935393 
1.9546477 
2 o2464576 
2o56l0418 


170.95484 

169 o 29360 

167 o 06122 

165.33070 

164.56022 

I64. 85080 

166.24802 

169.77120 

-165.60600 

-I64. 13332 

-163.50245 

-162.89356 

-162.55455 

-162.77557 

-163.13041 

-I63. 88891 

-165.12379 

-166.58995 

-168.18580 

-169.94772 

-171.82287 

-173.77752 

-175.78511 

-177.84262 

=-179.93537 
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247. 
248. 
249. 
250. 

251. 
252. 

253. 
254. 
255. 
256. 

257. 
258. 

259. 
260. 
261. 
262. 
263. 
264. 
265. 
266. 
267. 
268. 

269. 
270. 
271. 


2.8795196 
3.1896662 
3.4633445 
3.6690512 

3.8103315 
3.8769698 
3.9217889 
3.9790625 
4.1409762 
4.4609625 
4.9953985 
5.8880954 
7.2656523 
8.1811000 
.90562088 
.81687892 
.76044280 
.71118680 
.68455013 
.69116844 
.71238788 

.76123619 
.84674989 
.96413752 
1.1133331 


177.95661 
175.87351 
173.86999 
171.98936 
170.31615 
168.82943 
167.55483 
166.78406 
166.88418 
167.91345 
169.71161 
171.96057 
173.31975 
175.57659 
-163.42604 
-162.22822 

-161.48515 
-160.67153 
-160.23044 
-160.56392 
-161.12164 
-162.21030 
-163.85739 
-165.71847 
-167.65178 
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272. 
273. 
274  o 
275. 

276  a 

277. 
278. 

279. 
280. 
281. 
282. 
283. 
284. 
285. 
286. 
287. 
288. 
289. 
290. 
291  o 
292. 
293. 
294. 
295. 
296. 


1.3048910 
1.5428370 
1.8304559 
2.1665757 
2.5513942 
2.9773257 
3.4219972 
3.8637190 
4. 2571099 
4.5479984 
4.7205917 
4.7548153 
4.7096263 
4.6451493 
4.6905389 
4.9106693 
5o35246l0 

6.1322794 
7.3026670 
8.8962829 
.80854982 
.72256166 
.66983155 
.62448956 
.60207424 


-169.68291 

-171.75315 

-173.83722 

-175.92528 

-178.03303 

179.83285 

177.67307 

175.52600 

173-45832 

171.53646 

169.86121 

168.47200 

167.52009 

167. 17105 

167.78207 

169.35573 

171.62019 

174.05608 

176.00718 

177.29166 

-161.11039 

-159.74800 

-158.82945 
-157.74859 
-157.16038 
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297. 
298. 
299. 
300. 
301. 
302. 

303. 
302f. 
305. 
306. 

307. 
308. 
309. 
310. 

311. 
312. 

313. 
314. 
315. 
316. 
317. 
318. 
319. 
320. 
321. 


.61349231 
.64091441 
.69805652 
.79510355 
.92885725 
1.1011864 
1.3251106 
1.6076159 
1.9553803 
2.3711416 
2.8598067 
3.4176876 
4.0214149 
4.6391439 
5.1997527 
5.6082544 
5. 8143937 
5.7769834 
5. 5581712 
5.2677903 
5.1227914 
5.1854102 

5.4581457 
6.0624303 
6.9819H4 


-157.63016 
-158.42782 
-159.90782 

-162.02933 

-164.31382 

-166.57697 

-168.84837 

-171.07635 

-173.25640 

-175.40321 

-177.55745 

-179.75847 

177.96046 

175.63347 

173.35727 

171.24747 

169.42017 

167.94807 

167.00345 

166.69825 

167.60781 

169.68499 

172.70611 

175.88456 

178.84963 
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322  o 

323  o 

324  * 

325  c 
326o 
327a 
328. 

329  o 

330  o 

331. 
332  o 
333. 
334. 
335 . 
336, 
337. 
338. 
339. 
340  o 

3U. 

342  o 

343  • 

344  o 

345  o 
346o 


8,3882026 
.71704060 
063214025 
.58186488 

o 53962867 

o  52044645 

. 53503795 
.56648751 

.62899895 
• 73387204 
.87964441 
1.0703237 
1.3218132 
1.6447935 
2.0502029 
2.5459957 
3.1461223 
3.8584401 
4o6653009 
5.5299646 
6.3464236 
6.9395367 
7 .2020780 
7.O5OIOI4 
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NUMERICAL    RESULTS    OF    HAND   CALCULATIONS 

POINT  MAGNITUDE  (cm)      PHASE     (Degrees) 

1.  1.414  315° 

2.  9.86 
3-  10.1 

4.  1.414 

5.  7 

6.  14 

7.  3 

8.  11.4 
9c  6 

10.  3.14 

11.  7.07 
12  o  6.0                                      0 

13 .  4.0  0 

14.  4.47  60° 

15.  7.07  8.14° 

16.  11.05  5.2° 


24" 

0 

354 

45° 

0 

0 

0 

343.3° 

90° 

14° 

0 

8.14 
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NUMERICAL    RESULTS    OF    16    POINT    GRID 
COMPUTER    SOLUTION 


POINT 

AMPLITUDE    (cms.) 

PHASE    (DEGREES) 

1. 

3.47 

27.8 

2. 

3.43 

37.35 

3. 

3.16 

45.19 

4. 

3.16 

56.61 

5. 

4.60 

13.02 

6. 

4.94 

14.65 

7. 

4.95 

17.13 

8. 

3.74 

32.58 

9. 

4.18 

49.20 

10. 

5.72 

9.H 

11. 

6.37 

8.11 

12. 

8.28 

3.06 

13. 

3.04 

2*98 

14. 

5.54 

5.41 

15. 

6.49 

7.76 

16. 

6.98 

7.30 
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